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Abstract 

The energy losses of heavy multiplycharged ions at collisions with light atoms and polarization 
losses at moving through matter have been considered under circumstances the ion charge Z>1 and 
the relative colliding velocity « > 1, so that Z ~ v < c , where c - the light velocity (atomic units are 
used). In this region of parameters the Born approximation is not applicable. The simple formulas 
for effective stopping are obtained. The comparison with other theoretical results and experiments are 
given. 

1 Introduction 

Usually an energy loss of a relativistic charged particle in collisions with an atom is calcu- 
lated in the Born approximation which is applicable if Z/v <C 1 (here Z - the ion charge, 
v - relative velocity of collision, atomic units are used) (Berestetskii et al 1989). However, 
there were some resent experiments which use ions of such a large charge so the above 
mentioned condition was violated even at v < c ~ 137 (see for instance Kelbch et al 1986, 
Berg et al 1988, Berg et al 1992, Scheidenberger et al 1994). The approximations appli- 
cable at Z/v ~ 1 - Eikonal approximation and its modifications (McGuire 1982, Crothers 
and McCann 1983), the Method of Sudden Perturbations (Eichler 1977, Salop and Eich- 
ler 1979, Yudin 1981, Yudin 1991), as well as the Method of the Classical Trajectories 
(Olson 1988) demand considerable numerical calculations even in nonrelativistic range of 
collision velocities. Relativistic collision velocities make the the calculations even more 
complicated (Yudin 1991). 

In the present paper we consider the energy loss of heavy relativistic highcharged ions 
in collisions with light (nonrelativistic) atoms and polarizations losses at moving through 
the matter under conditions Z ~ t> < c, Z^>l,f^>lby following a simple approach 
(Matveev 1991) and its relativistic modification (Matveev and Musakhanov 1994). We are 
going to obtain simple formulas describing an effective stopping and compare it's results 
with the results of other theoretical methods and experiments. 



2 Collisions with separated atom 



Accordingly to (Landau and Lifshitz 1974) average losses of energy in collisions are char- 
acterized by the effective stopping 



K 



Here e n and eo are the energies of excited |n > and bound |0 > atom's states, a n is the 
cross section for excitation of state \n >. First, we consider, for simplicity, the collision 
of a relativistic highcharged ion with a hydrogen atom. Accordingly to (Matveev 1991, 
Matveev and Musakhanov 1994) the whole interval < b < oo of all possible values of 
impact parameter b may be divided into three ranges: 

A) < b < bt ; B) b x < b < b ; C) b < b < oo , (2) 

which correspondent to small, middle and large values of an impact parameter. The 
border amounts b\, bo are (Matveev and Musakhanov 1994) 



6 ~ v/sjl-F- (3 = v/c. 



We are going to calculate k in each range from (Q) and obtain total effective stopping 
by summing those values . The exact meaning of bo and b\ are unimportant for us, since 
dependence k upon sewing parameters b\, bo appears to be logarithmic, which lead to the 
correct sewing of contributions from each region (^) and vanishing any dependence k on 
b\ and bo in the final result. 

The range of small impact parameters is the range of large momentum transfers, thus 
we needn't take into account the interaction of an electron with an atom nuclei, and 
therefore we can consider the scattering of an ion on an electron having rested before 
collision (Berestetskii et al 1989 §82). In this case an effective stopping may be expressed 
through using a(e) which is a cross section for collision with energy transfer e: 

tax 

ea{e)de . (3) 
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Here we can't, unlike (Berestetskii et al 1989), use a(e) in the Born approximation, because 
a condition Z/v <C 1 is invalid. To go on, we assume colliding ion to be an infinitely heavy 
particle, and consequently it doesn't change its moving. This case the cross section for 
scattering of an ion moving at fixed velocity on the free and motionless electron may 
be obtained by transformation to the system where ion is motionless. We mark 9 a 
respective angle of scattering. Following to (Berestetskii et al 1989) we can maintain an 
energy transfer to be a function e = e(#), except the case of a superhigh energy. That 
function has a form 



The values e max and e min from (|) are reached at 6 = it and 6 = 9 min respectively. As a 
result we can rewrite (|3|) in a form 

v 2 Je min o R {6) 2 

Here a (6) is the cross section for the scattering of electron on motionless ion with a 
charge Z (at electron velocity v), which is obtained (Ahieser and Berestetskii 1969, Mott 
and Massey 1965) from exact solution of a scattering problem for the Dirac equation; or 
is the Rutherford cross section: 

Z 2 (l-/3 2 ) 

aR c A p 4 (i-cosey ' K ) 

The ratio a{6)/an{6) — > 1 at 9— >0 (Ahieser and Berestetskii 1969, Doggett and Spenser 
1956), and hence in order to define Q m in we can use the quassiclassical connection (Landau 
and Lifshitz 1973) between a scattering angle and an impact parameter (see also the 
obvious qualitative picture of collision suggested in (Matveev 1991)) 

9 mn = ^-Jl-P 2 ■ (7) 
v 2 bi v 

One can easy see that the integral fl5|) depends on angle 9 min in logarithmic way (at small 
9 m in) , thus the formula @ may be written in a form: 

« = 47r^ln H - -, (8) 



v 



Zy/1 - P 2 a(Z, v) 



The function a(Z,v) are defined by us from comparing (|S[) with the numerical calculation 
results by formula (El) with the ratio ct/cxr from (Doggett and Spenser 1956). As a result 
we can approximate function a(Z,v) with a following formula: 

a{Z, v) = (-0.23016 • a - (1.00832 • a - 0.32388)/? 2 + l) 2 , (9) 

here a = Z/c. In the tables number 1 and number 2 the effective stopping calculated (at 
bi = 1) by formula (|]) (table 1) and by formula (§) with substituting the function a(Z,v) 
from (|9|) (table 2) are brought. In these tables the data of 1st column are ion's energies (in 
Mev/nucleon) correspondent to the relative velocities from (Doggett and Spenser 1956); 
the data of next columns are k values (in atomic units) for the ions with a charge: 6, 13, 
29, 50, 82, 92 respectively. One can see from the tables that the supposed approximation 
is enough good, at any rate in the limits of v and Z which the data of (Doggett and 
Spenser 1956) are brought for. Besides we have to emphasize that function a(Z, v) — > 1 
at (3 — > (nonrelativistic limit) and a — > 0, which conforms to the fact that o jo^ — > 1 at 
(3 -»• 0, a -> 0. 

The collisions with the middle impact parameter b\ < b < 6 make the main con- 
tribution to the cross sections for the inelastic process (Matveev 1991, Matveev and 
Musakhanov 1994). The energy transfer in that collision is e < 1 ~ ionization potential 
of an atom, and therefore the contribution to the effective stopping from that collision 
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Tabic 1: Effective stoping (at. units) - result of numerical integration by formula (5) with using 
from (Doggett and Spenser 1956) 
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Table 2: Effective stoping (at. units), by formulas (8) with using a funstion a(Z,v) defined from 



can't be accounted by using a perturbation theory (Eichler 1977, Yudin 1981, Matveev 
1991, Matveev and Musakhanov 1994). We must note else, that an atom electron is non- 
relativistic as before so after collision with an impact parameter from range b\ < b < b Q 
(Matveev and Musakhanov 1994). The contribution of a collision with such impact pa- 
rameter to an effective stopping may be easy obtained from formula ([!]) by substituting 
to it the cross section for the inelastic process (Matveev and Musakhanov 1994) 

o n = / 2nb \<n\ exp(-iqf) | >| db , (10) 

J b± 

Here q = 2Zb/(vb 2 ) and b is impact parameter vector, r is a coordinate of atomic electron. 
On repeating the calculations of (Landau and Lifshitz 1974) (but our case is simpler, 
because the upper limit in integration (|T^) doesn't depend on the final state of an atom) 
we have for an effective stopping: 

Z 2 qi 

K = 5^(^n Co)^n = 471 — In , (11) 

n v 2 q 

here q = 2ZJ (vb ); qi = 2Zj [yb\) . 

We need to account the contribution from collisions with an impact parameter belong- 
ing to the range bo < b < oo. In this range the interaction between an ion and an atom 
may be accounted by using a perturbation theory (Eichler 1977, Matveev 1991, Matveev 
and Musakhanov 1994). Amplitude for the transition of an atom from a state |0 > to a 
state | n > may be obtained following to (Moiseiwitsch 1985) 
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here Q n = e n — e ; £ = VL n byJl — (3 2 /v ; K (£), i^i(£) - the McDonald functions, 
ron =< 0|f|n >. 

The cross section correspondent to (|P2"D is 



a n {b < b < oo) = J d 2 b \ A 0n 



|2 



The integrating in this formula is made in limits: the angle of vector b changes from to 
2tt and 6o < b < oo. In a result the cross section has a form []: 

Z 2 ( 4v 2 \ 

° n = An V 1 Xon |2 ( ln v 2 b 2 nui - p 2 ) ~ ' (13) 

Here rj = e B = 1.781 (B=0.5772 - the Euler constant), XQ n =< n\x\0 >. The contribution 
to a effective stopping from collisions with large impact parameters may be calculated 
through substituting expression (|13D to formula ([!]) 



1 To say rigorously, formula (|l^) was obtained from supposition that § = O n 6^/ 1 — ffl Jv <C 1 (f2 n ~ 1), then b <g 



1 — /3 2 ~ bo, so following sewing may be performed exactly at such b. 



Here we have introduced "average atom energy"- I (Berestetskii et al 1989 §82), which is 

j n j _ ^"( 6w ~ £ °) I X ° n I 2 m ( e " ~ e o) 



The total effective stopping of a relativistic highcharged ion on a hydrogen atom is 
obtained on summing the results of formulae (|l4]),(|l"l~l) and @: 

Z 2 { 2v 3 1 

" =4 v| Vz(i-w^) -^)- (16) 

We would like to bring k value calculated in the Born approximation from (Berestetskii 
et al 1989) 

2v 2 

v 2 { 11 /(l-/? 2 

It should be noticed that nonrelativistic limit for ionization losses flTSl ) (in accounting 
of a(Z, v) — > 1 at f3 — > and a — ► 0) 



4 ^ In— — -/? 2 . (17) 



Z 2 2v 3 

k = Atc— In-— (18) 

has a form which is similar to wellknown Bohr formula (Bohr 1913) obtained from classical 
representations. Besides we have to notice that our result (|16|) obtained on the basis of 
approach (Matveev 1991, Matveev and Musakhanov 1994), which is valid at Z ~ v ^> 1 
only, doesn't have a way to pass (so does the Bohr formula) to the Born approximation 
which is valid at v/Z -C 1. 

As followed from our way of obtaining formula (0), in order to generalize it as well as 
the Born losses (|T7|) to the case of collisions with multielectron atoms (whose electrons 
have velocities v a <C v, v is ion's velocity) we need to multiply the right part of formula 
fll6|) on Z a [Z a is the number of atom electrons) and change value I for average atom 
potential l a defined, as before, from (|i"5l), but with values e n , e and | x 0n \ 2 calculated for 
complex atom. 



3 Energy losses in matter 

Let's consider an energy loss of relativistic highcharged ion moving through matter. This 
losses are the sum of macroscopic (polarization) losses and losses at collision with sepa- 
rated atoms. The ion velocity hold to be much bigger than typical velocities of atomic 
electrons or at least the majority of them. Accordingly to (Fermi 1940) polarization losses 
of a charged particle moving in the matter are defined by the flow of the energy of elec- 
tromagnetic field of a particle through the cylinder with radius b' Q built around a particle 
trajectory. The effective stopping is obtained by dividing of the flow on a particle velocity: 

Z 2 h' r+°° ( 1 \ 

k = —±\ duK Q {b'SK x {b' o e) [-TT - P 2 ) , (19) 



here £ 2 = oj 2 (v 2 — c 2 s(uj)), e(oj) - dielectrical penetration coefficient. At enough small 
cylinder radius b' , i.e. 

I b' £ |« 1 (20) 

we have 

IZ 2 r+°° f 1 1 2 

«= — / udul-j-^-p 2 }]*—, (21) 



where 77=1.781 (like (|13D ) . Other side, accordingly to (Landau and Lifshitz 1982) energy 
loss may be counted as a work of electromagnetic field per unit of way 

K = / / uduqdq - \ >> . 22 



Formula (21) will coincide with ( P2j) if condition (|20| ) is valid and go — 2/ (t7& ): surely, 

2 1 , / 2 2 \ rio qdq 

There are two distinguished cases to be considered separately: a) v 2 < c 2 /sq (eq = e(0)) 
- dielectrical penetration coefficient in static field) and b) v 2 > c 2 /eo- In first case (Landau 
and Lifshitz 1982): 



K = An- 
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q = 2/(r76 ) . (24) 



u^/l^W 2 

Here uo is a average amount of frequency of atom's electrons motion: 

f£°ujr)"(u}) In uj du 

hxu= ' \ , 

Jo ur]"(Lj)dLJ 

here 77" = Im{e^ 1 {uj)). In second case ((v 2 > c 2 /eq))\ 1) if the particle energy is not 
rather large (motion energy is either less than or order of a rest energy of ion) we can use 
the formula (|23|); 2) In ultrarelativistic case (Landau and Lifshitz 1982) 

^2^1n^. (25) 
c 2 AixN K ' 

Now we must sum the macroscopic losses and the energy losses on separated atoms. 
In order to do it rewrite the condition (|20T) in a form: 



b' < vj I uyjl- (3 2 s \< v/^/l-(3 2 ~ b 



where, for estimations, was assumed that u ~ ou a ~ 1 - typical atomic frequency. The 
condition of macroscopic method's being applicable is b' ^> bi ~ 1 - typical atomic size. 
Thus the lowest border b' is between b\ and b : 

h < b' < 6 (26) 

On comparing this condition with (^j we conclude that in order to obtain the total 
energy losses of relativistic highcharged ion moving through the matter we have to sum 
the contribution from ranges A and B from (0) (in range B the upper limit is equal to 
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Table 3: Experimental and theoretical meanings of energy losses (in M ev j {rag / cm 2 )) for various combi- 
nations target-ion. 
(*)-Scheidenberger et al 1994 

b' Q ) and polarization loss. The sum of contribution from ranges A and R according to (|8|) 
and (|Tl"D is equal 



Z 2 N , 
47r — — In 

V 2 



v 2 U 



Za{Z,v)^l- 



(27) 



We changed the number of atom electrons - Z a for the number of electron in the unit 
of volume - N (as was hold in Landau and Lifsitz 1982). Ry summing ( |2"3| ) and (^) we 
obtain the total energy loss of relativistic highcharged ion moving through matter for the 

case v 2 < c 2 /eo 
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Zt](1 - (3 2 )a(Z, v)Q 

As was said this formula is used often for the case v 2 > c 2 /eq and for the particle whose 
velocity is not too high. It should be noticed that (|22"D is different from (|18|) describing 
the energy losses on separated atoms with changing average potential I for Q only (there 
is the same situation in Landau and Lifshitz 1982). In ultrarelativistic case on doing the 
same way as before we obtain total effective stopping in a form: 
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(29) 



Z 2 r] 2 (l- (3 2 ) 2 a 2 (Z,v)4irN 

One can see that in this case accounting the polarization loss leads to slower growing an 
effective stopping power with increasing a collision velocity, in comparing with case of 
losses on separated atoms (0). One may rewrite the formula (|29|) introducing the Fermi 
corrections on density effect: 
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(30) 



Z V (l-f3 2 )a(Z,v)I a 

here the meaning of average ionization potential I a and the Fermi corrections 5/2 are 
brought from (Sternheimer et al 1984, Inokuti and Smith 1982). 

In table 3 the theoretical and experimental meaning of effective stopping 
(in Mev / {mg / cm 2 )) for various ions (we limited ourselves with Kr and Xe which have 



enough large charge) and targets from (Scheidenberger et al 1994), as well as our results 
are brought: the column 1 - ion, it's energy (in Mev/nucleon) and amount (3 = v/c; the 
column 2 - the target's nature; the column 3 - the results of calculations (Scheidenberger 
et al 1994) by the Bete formula with the Fermi corrections; the column 4 - the results of 
calculating (Scheidenberger et al 1994) by the Bete formula with accounting Fermi correc- 
tions, Mott corrections (Mott 1929) and Bloch corrections (Bloch 1933); column 5 - are 
our results; column 6 - the experimental meaning for effective stopping (Scheidenberger 
et al 1994). One can see that our results are in rather good accordance with experiment. 

4 Conclusion 

The simple approach suggested in the papers (Matveev 1991, Matveev and Musakhanov 
1994) enable us to estimate an effective stopping of relativistic highcharged ion in col- 
lisions with separated atoms and at motion in matter for many practically important 
cases, since the formulas obtained in this paper enable us to use the wellknown method 
of phenomenological corrections usually used in applied calculation. The region of appli- 
cability of our formulas Z ~ v < c doesn't permit to make direct transition to the Born 
approximation (Z/v <C 1). But that fact is not too significant limitation, because the 
region of applicability of the Born approximation for the ion with enough large charge 
(for example Z = 92) is unreachable even at v ~ c. 

The authors acknowledge Dr. Cristoph Scheidenberger, GSI, Darmstadt, Germany for 
attention to our work and information about his results. 
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